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Continuing previous work reported in an earlier paper [L.M. Burko, A.I. Harte, and E. Poisson, 
Phys. Rev. D 65, 124006 (2002)] we calculate the self-force acting on a point scalar charge in 
a wide class of cosmological spacetimes. The self-force produces two types of effect. The first is 
a time-changing inertial mass, and this is calculated exactly for a particle at rest relative to the 
cosmological fluid. We show that for certain cosmological models, the mass decreases and then 
increases back to its original value. For all other models except de Sitter spacetime, the mass is 
restored only to a fraction of its original value. For de Sitter spacetime the mass steadily decreases. 
The second effect is a deviation relative to geodesic motion, and we calculate this for a charge that 
moves slowly relative to the dust in a matter-dominated cosmology. We show that the net effect of 
the self-force is to push on the particle. We show that this is not an artifact of the scalar theory: 
The electromagnetic self-force acting on an electrically charged particle also pushes on the particle. 
The paper concludes with a demonstration that the pushing effect can also occur in the context of 
slow-motion electrodynamics in flat spacetime. 



I. INTRODUCTION 

The prospect of a near-future measurement of gravi- 
tational waves produced by the capture of a solar-mass 
compact object by a massive black hole (one of the iden- 
tified sources for the Laser Interferometer Space Antenna 
0) is motivating an intense effort to calculate the mo- 
tion of a compact body in a small-mass-ratio approxima- 
tion that goes beyond the test-mass description. In this 
approximation the gravitational influence of the smaller 
body is taken into account, but the motion is still de- 
scribed as taking place in the background spacetime of 
the larger body. While the motion is geodesic in the test- 
mass description, it is accelerated in the small-mass-ratio 
approximation, and the small body is said to move un- 
der the influence of its gravitational self-force. The work 
that has gone into formulating and evaluating thegravi- 
tational self- force is reviewed in detail in Refs. @J3| , and 
a shorter, less technical account is given in Ref. [J. 

The calculation of the gravitational self-force involves 
a number of technical and conceptual difficulties, and 
concrete results pertaining to situations of astrophysical 
interest have not yet been obtained. To shed light on the 
physics of self-forces in curved spacetime it is helpful to 
consider toy problems. In this paper we continue an ef- 
fort begun in an earlier paper (Ref. @, hereafter cited as 
"paper I" ) to study the self- force acting on a point scalar 
charge inserted in a cosmological spacetime. The self- 
force acting on a scalar charge displays many similarities 
with the gravitational self-force (including a dependence 
on the particle's past history), but it is technically much 
easier to obtain (fewer field components are involved). 
And the cosmological models considered here are suffi- 
ciently simple that the scalar self-force problem can be 
solved in a countably infinite number of examples. Be- 
yond these technical advantages, the physics of a scalar 
point charge is surprisingly interesting: Because a scalar 
field can radiate in monopole waves, the inertial mass of 
a scalar charge is not necessarily a constant of the mo- 



tion. This phenomenon of mass radiation is one of our 
main topics in this paper. We also calculate the motion 
of a scalar charge when it is subjected to its self-force, 
and compare this with the motion of an electric charge. 

The physics of scalar charges in curved spacetime was 
first investigated by Quinn [a; his theory is briefly re- 
viewed in Sec. [HJ In paper I [5| Quinn's theory was ap- 
plied to a scalar particle at rest in a spatially-flat cosmo- 
logical spacetime with metric 

ds 2 = a 2 (r])(~drf + dx 2 + dy 2 + dz 2 ), (1.1) 

where rj is conformal time and (x, y, z) are comoving spa- 
tial coordinates. In paper I the scale factor was restricted 
to the forms a(rj) — Cr/ 2 (a matter-dominated cosmol- 
ogy) and a(r/) = C/\rj\ (de Sitter spacetime); C is a scal- 
ing constant. For the matter-dominated cosmology it 
was shown that the scalar charge initially loses a fraction 
(possibly all) of its mass, but that this radiated mass 
eventually returns to the particle. For de Sitter space- 
time it was shown that the particle radiates all of its 
mass within a finite proper time. 

In this paper we generalize the work reported in paper 
I 5] on several fronts. We consider a much larger class 
of cosmological models, those for which a(rj) = C\rj\ a 
with a restricted to (positive and negative) integers. For 
this class of models we obtain an exact expression for 
771(77), the particle's mass function. We show that for all 
values of a except a = — 1 (de Sitter spacetime), the 
mass first decreases and then increases back to a fraction 
of its original value; when a is positive this fraction is in 
fact equal to unity. 

We also calculate the self-force acting on a scalar 
charge that moves relative to the cosmological fluid. Here 
our treatment is no longer as general: it is restricted by a 
slow-motion assumption and to a matter-dominated cos- 
mology. We show that the net effect of the self-force is 
to push on the particle: When a charged particle and an 
uncharged particle start with the same initial velocity, 
the charged particle subsequently moves faster than the 
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uncharged particle; ultimately both particles are brought 
to rest by the cosmological expansion. While this con- 
clusion may appear surprising, we show that it is not an 
artifact of the scalar theory: the same conclusion applies 
to an electrically charged particle. 

We begin in Sec. [H] with a description of our cosmolog- 
ical spacetimes and a review of Quinn's theory of point 
scalar charges. The time-changing mass and the motion 
of the particle can both be determined once the retarded 
Green's function for the scalar wave equation has been 
computed; this computation is carried out in Sec. lIIII In 
Sec. IIVI we examine the inertial mass of a scalar charge 
at rest in our cosmological spacetimes, and establish the 
behaviors that were described previously. In Sec. we 
calculate the motion of a slowly-moving scalar charge in 
a matter-dominated cosmology, and show that the net 
effect of the self-force is to push on the particle. In 
Sec. IVII we carry out a similar calculation for an electri- 
cally charged particle, but we generalize the cosmological 
models to the class described by a (77) = Crf* for any (real 
or integer, but positive) a; we show that the electromag- 
netic self-force also pushes on the particle. In our final 
section fVll we reflect on this perhaps surprising property 
of the self-force, that it accelerates the particle instead of 
slowing it down (relative to the motion of an uncharged 
particle); this goes against a naive expectation that the 
self-force should always give rise to radiative damping. 
In Sec. IVII we show that this property is not an artifact 
of curved spacetime: It can be manifested also in the 
simplest context of slow-motion electrodynamics in flat 
spacetime. 



II. EQUATIONS OF MOTION OF A SCALAR 
CHARGE IN A COSMOLOGICAL SPACETIME 



The mass density of the cosmological fluid is 
3a 2 1 



8ttC 2 M 2a + 2 



and its pressure is given by 



3a 



(2.3) 



(2.4) 



Imposing the dominant energy condition (see, for exam- 
ple, Sec. 9.2 of Ref. 3 or Sec. 2.1 of Ref. ||) restricts a 
to the interval a > h if it is positive, and a < —1 if it 
is negative. For a = 0we have flat spacetime, for a = 1 
a radiation-dominated cosmology (p = -|p), for a = 2 a 
matter-dominated cosmology (p — 0), and for a = — 1 
we have de Sitter spacetime (p = —p). 

The scalar field <!> produced by the scalar charge q 
obeys the wave equation 



(2.5) 



in which V™ indicates covariant differentiation and 



n(x) = q / 8<t(x, z(t)) dr 



(2.6) 



is the scalar charge density; here 8i{x,x') is a scalar- 
ized, four-dimensional Dirac distribution that satisfies 
J 5i(x, x')y/—g d A x — 1 when x' is within the domain of 
integration. The physical solution to Eq. H2.5J) is written 
as 



$(x)=q / G(x,z(r))dT, 



(2.7) 



in terms of a retarded Green's function G(x 1 x') which is 
a solution to 



A particle of mass m and scalar charge q moves in a 
cosmological spacetime with metric 



ds 2 = a 2 (f])(-drf + dx 2 + dy 2 + dz 2 ) , 



(2.1) 



where r\ is conformal time and x = (x,y, z) are comov- 
ing spatial coordinates; we use x = (r\. x) to label a 
spacetime event, and boldfaced symbols represent three- 
dimensional vectors in flat space. The cosmology is re- 
stricted to be spatially flat, and the scale factor is re- 
stricted to the form 



g a ^V a VpG(x,x') = -4ttS 4 (x,x'). 



(2.8) 



When viewed as a function of x' , the Green's function 
has support both on and within the past light cone of 
the spacetime event x; the contribution from the light 
cone is singular, but the contribution from the interior is 
smooth. 

The field §(x) is singular on 7 but its action on the 
scalar charge is nevertheless well defined. As shown by 
Quinn (and reviewed in Refs. @ and 0), the equa- 
tions of motion of a freely moving scalar charge are 



a( V ) = C\n\ 



(2.2) 



where C is a scaling constant and a is a constant (positive 
or negative) exponent. When a > we take 77 to be 
positive and increasing from to 00. When a < we 
take r\ to be negative and increasing from —00 to 0. In all 
cases a(rj) increases with increasing 77. The particle moves 
on a world line 7 described by the parametric relations 
z M (r), with r denoting proper time; its velocity vector is 
u» = dz^/dr. 



m^— = q 2 {5^ + u^u y ) 



and 



dm 
~d7 



dr 



1 

12^ 
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V V G(z{T) lZ {T'))dT' 



(2.9) 



<" ' \»G(z{T),z{T'))dT' 
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(2.10) 
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The left-hand side of Eq. (|2.9(l is the product of m(r) 
with the particle's covariant acceleration vector, and the 
right-hand side is the scalar self-force. This involves the 
spacetime's Ricci tensor and the particle's velocity vector 
evaluated at the current position z(r) on the world line. 
It also involves the gradient of the retarded Green's func- 
tion evaluated at z(t) and a prior position z(t') — the 
Green's function is differentiated with respect to x before 
making the identification x — z(t). Finally, the self- force 
involves an integration over the past portion of the par- 
ticle's world line, from the big bang at t' = to almost 
t' = r; the integration is cut short by e = + to avoid the 
singular behavior of the Green's function at coincidence. 
The right-hand side of Eq. H2.10|) gives the rate at which 
the particle's inertial mass changes with proper time r; 
it involves the same objects as Eq. (|2.9|) . together with 
1Z = the spacetime's Ricci scalar evaluated at z(t). 

To solve Eqs. (|2.5|l and (|2.8|l in the cosmological space- 
times it is useful to introduce the reduced field 

^(r?,x) = a(r))$(x), (2.11) 

the reduced density fj,*(ri,x) — a 3 (rj) /i(ai) , and the re- 
duced Green's function 

5(77, x; r,', x') = a(rj)a(r)')G(x, x'). (2.12) 

These satisfy 

+ * + --*><• < 2 ' 13 > 

and 

( d 2 o a(a — 1)\ _. ,. „ . ., 

{-g^ + V + ^T^y = - ^ X - X )' 

(2.14) 

In these equations V 2 stands for the Laplacian opera- 
tor in three-dimensional flat space, and £3(2; — x') = 
8(x—x')8(y—y')5(z — z') is an ordinary three-dimensional 
Dirac distribution. The reduced density obtained from 
Eq. lE3 is 

fx*(r,,x) = %(x-z(r,)), (2.15) 

7 

where 7 = a(ri)dr]/dT = (1 — v ■ t?)~ 1//2 , with v(rf) = 
dz/d-q denoting the particle's coordinate velocity vector 
and v ■ v = (v x ) 2 + (v y ) 2 + (v z ) 2 . The retarded solution 
to Eq. l(2~T3) is 

#7,a:) = J g(^x;r } ',x')^(r ] ',x')di 1 l d 3 x', (2.16) 

and this statement is equivalent to Eq. I|2.7|) . 

III. REDUCED GREEN'S FUNCTION 

The difficulty in evaluating Eqs. I|2.9|l and l|2.10|l re- 
sides mostly in the computation of the retarded Green's 



function. In paper I |5j] the Green's function was ob- 
tained for two specific cosmological models, those for 
which a — 2 (matter dominated) and a = — 1 (de Sit- 
ter). In this section we generalize this calculation and 
compute the Green's function for arbitrary, but integer, 
values of a. 

For convenience we introduce a positive integer I such 
that 

a = — I when a < (3.1) 

and 

a = I + 1 when a > 0. (3.2) 

We notice that with this substitution, a(a — 1) = 1(1 + 
1). We notice further that a given value of I produces 
two distinct cosmological models, one with a\ = — I < 
and another with o?2 = I + 1 > 0; these models are 
"equivalent" in the sense that ai(ai — 1) = 0:2 (c*2 — 
1) = 1(1 + 1), so that they give rise to the same reduced 
Green's function. In this sense, flat spacetime (a = 0) 
is equivalent to a radiation-dominated cosmology (a — 
1), and dc Sitter spacetime (a — —1) is equivalent to a 
matter-dominated cosmology (a = 2). 

It was explained in some detail in paper I p| that for 
the cosmological spacetimes considered in this paper, the 
general structure of the reduced Green's function is 

g(r),x;r]',x') = 

+ %-r/ -R)V(v,x;ri',x'), (3.3) 

where R = \x — x'\ = \J (x — x') 2 + (y — y') 2 + (z — z') 2 
and V(j], x; r/, x') is a smooth two-point function that 
satisfies the homogeneous form of Eq. I|2.14|l . The first 
term on the right of Eq. (|3.3|) is the singular contribution 
from the light cone; the second term is the contribution 
from the cone's interior. In flat spacetime V — 0, and the 
support of the Green's function is entirely on the light 
cone. Because a radiation-dominated universe (a = 1) is 
equivalent to flat spacetime (in the sense defined in the 
previous paragraph), we have that V = in this case as 
well. 1 In the following we will exclude / — from our 
considerations, and we will compute the reduced Green's 
function for r\ — rf > R; in this domain we have g = V. 

A Fourier-transform method to solve Eq. (|2.14|) was 
described in Sec. IV B of paper I :5|, and the reader is 
referred to this paper for details. The core of the method 
consists of integrating the ordinary differential equation 

( d 2 2 Z(Z + 1)\ . . 

[d^ +k - - V>=° (3 - 4) 



This conclusion follows also from the fact that 72. = for a 
radiation-dominated cosmology, which implies that the scalar 
wave equation is then conformally invariant. Because the space- 
time is conformally related to flat spacetime, the cosmological 
and Minkowski solutions are related simply by a conformal fac- 
tor. 
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for the function (7(77, 77'; k), subjected to the boundary 
conditions 

.9(77 = 77'; fc) = 0, ^(77 = rj';k) = 4?r. (3.5) 

The retarded Green's function is then obtained by invert- 
ing the Fourier transform, which leads to 

1 f°° 

V = —r— / g(r),r] ;k)ksm(kR)dk 
2n 2 R J 

1 d f°° 

= - 2 TT 2 RdR J 9(v,V,k)cos(kR)dk (3.6) 

when 77 — 77' > R. The second representation will be more 
useful for our purposes. 

The general solution to Eq. 13. 411 is g = Akrjji(k\rj\) + 
Bkr)ni(k\r)\), which is expressed in terms of spherical 

I 



Bessel functions. The boundary conditions of Eqs. i|3.5|) 
determine the constants A and B, and their form can be 
simplified by invoking the Wronskian-type relations that 
are satisfied by the Bessel functions (see, for example, 
Sec. 10.1.6 of Ref. @). We find A = -47r|77'|n ; (fc|7/|), 
B = 4ir\r]'\ji{kW\), and 



^Awkvlv'llMklv'l^ikl^-miklv'lMk^)]. (3.7) 



After substituting this result into Eq. (|3.6|l and us- 
ing the identities (see, for example, Sees. 9.1.2, 10.1.1, 
and 10.1.12 of Ref. @) ji(z) = {ir/2z) 1 / 2 J l+ i(z), 

ni(z) = (-1Y +1 (tt/2z) 1 ^J_ 1 _i(z), and cos(z) = 

{itz/2) 1 /' 2 J_ 1 (z), we obtain 



T/ = sign(«)(-l) f 



2R8R 



VWR / V^r^'i+i(*l»7 / l)-'_i-i(fcl»?l) - ^-i_i(*l»7 , ])J'i + i(*l»7l)V-i(*-H)d*: s (3-8) 



in which all functions are Bessel functions of the first kind. The expression also involves sign(a) = sign(?7) = r)/\r}\; 
recall that 77 and 77' are both negative when a < 0. 

The integrals on the right-hand side of Eq. (|3.8(l can be evaluated in closed form, thanks to Eq. (6.578.1) of Ref. [To| . 
In the next paragraph we shall describe the calculation of V for a > only; a very similar calculation that would 
return V for a < will not be presented. 

The first integral on the right-hand side of Eq. (|3.8|) is proportional to l/r(— I) and it vanishes when I is an integer. 
The second integral is equal to 



(v/v 



r(-j + i)r(j + i) V wtfR 



:F 4 



-l,±;-l + l±;(v'/v) 2 ,(R/v) 2 } 



and it is expressed in terms of Appell's hypergeometric 
function of two variables (see, for example, Sec. 9.18 of 
Ref. 01)- Because its first argument is a negative in- 
teger, the Appell function is in fact a terminating poly- 
nomial in the two variables (rj'/r/) 2 and (R/rj) 2 . After 
cleaning up the Gamma functions we arrive at 

(2/ -1)11(77/77')' 



V = 



R 



x [-1, |; -I + |, |; (77777)2, (R/ V ) 2 ] (3.9) 



(20!! 
d_ 

OR 

for the reduced Green's function. There remains the task 
of evaluating the derivative with respect to R. We use 
the identity 

—F 4 (a,(3;-f,j';x,y) = —F 4 (a+l,(3+l;j,j'+l;x,y), 
oy 7 

which is easily established from the series representation 
of the Appell function. Our final result is 



V(r],x;i]',x') 



(2i-l)!! (77/77') 
(2Z-2)!! T] 2 



x F 4 [-1 + 1, |; -I + i, |; (77'/77) 2 , (R/ V ) 2 ] , (3.10) 

where R = \x — x'\. The complete reduced Green's func- 
tion is then given by Eq. (|3.3() . and G(x,x') can be ob- 
tained from Eq. I|2.12|l . 

The calculation for a < returns a different ex- 
pression for V, but this can be obtained directly from 
Eq. H3.10fl by performing the exchange r\ <-> 77'. These 
expressions are in fact identical, because the two-point 
function satisfies the reciprocity relation V[r\\ x'; 77, x) = 
V{r],x\r[,x'). This property is derived, for example, in 
Refs. HQ. We conclude that Eq. j3~TU)) gives the correct 
expression for V for both positive and negative values of 
o. 

The Appell function reduces to an ordinary hyper- 
geometric function when the second argument is zero: 
^4(0;, P; 7, 7'; x, 0) = F(a, P;~/;x). This implies that 
Eq. H3.10|l simplifies to 



V(r],x;r}',x) 



r\l 



(2Z-1)!! (77/77') 
(2Z — 2)!! 7 ? 2 
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xF[-l + -I + i;(7 ? 77 7 ) 2 ] (3.11) 

when the spatial points x and a;' coincide. When the 
spacetime events x and x' fully coincide we have 



V(r],x;r),x) = 



1(1 + 1) 
2r/ 2 : 



(3.12) 



which is recovered after using identity (15.1.20) of 
Ref. This result, when combined with Eqs. (|2.12|l 
and 12.2fl . implies that the full retarded Green's function 
has the coincidence limit 

G( V , x; V -e,x) = a ^ 1 ^ a+2 = ^ft, (3.13) 

where 1Z is the spacetime's Ricci tensor at x = (77, x). 
The presence of e = + in Eq. (|3.13(l serves to remind 
us that the coincidence limit excludes the singular, light- 
cone contribution to the Green's function; this e plays 
the same role in Eqs. I|2.9I) and l|2.10|l . 

The general expression of Eq. (|3.10() reduces to simple 
forms when I is a small integer. We list a few in Table 
I. The result for I = 1 was first derived in paper I 
The general expression of Eq. i|3.10fl , and the additional 
special cases listed in the Table, are new. 



charge was created at t = To with an initial mass to(to). 
As a final step we change the variable of integration from 
t' to 77', using the differential relation dr' = a(rj')dri' 
appropriate for a particle at rest, and we re-introduce 
the reduced Green's function of Eq. I|2.12|l . Taking into 
account Eq. (|3.3|) we arrive at a mass function given by 



771(77) = m o — 



a(rf) 



V(v,0;r ) ',Q)dr 1 ', 



(4.2) 



'/(> 



where mo = 772(770) is the particle's mass at the moment 
of creation. 

A more concrete expression for 777(77) is obtained by 
substituting Eq. I|3.11[> into Eq. 14.21) . The resulting ex- 
pression simplifies when we change the variable of inte- 
gration from 77' to /i = (t)' /f]) 2 , which we recognize as the 
argument of the hypergeometric function. Simple manip- 
ulations yield 

777(77) 1 (2f- 1)!! +1 . 

= 1 ~ 7v — ^TjCwoA?) sign a 

mo (21 — 2)!! 

x / M -< z+1 >/ 2 F(-Z + l,§;-/ + ±;//)d//, (4.3) 



where 



IV. MASS CHANGE OF A CHARGE AT REST 

In this section we calculate the time-changing inertial 
mass of a scalar charge that is comoving with the cosmo- 
logical fluid. Such a particle moves on a geodesic of the 
spacetime, and as we shall see in Sec. E] the right-hand 
side of Eq. (|2.9|) vanishes in this case. The right-hand 
side of Eq. p.lOjl does not, however, and we shall find 
that the particle's mass changes with conformal time 77. 
Without loss of generality we place the scalar charge at 
the spatial position x = 0. 

We begin by deriving a convenient expression for the 
mass function. Equation (|2.10|l can be written as 



dm 
~dr~ 



1 

-2 U 



T ~ e d 

-j-G(t,t') dr' 
dr 



or as 

dm 
~dr~ = 



-U-G(T,T-e) 



G{T,T>)dT>. 



Taking into account Eq. Ij3.13|l we find that the first 
term vanishes, and that the equation can be directly in- 
tegrated. This gives 



m,(r) = m(r ) 



G^T^dT'. (4.1) 



We have truncated the domain of integration from t' = 
to t' = To > to avoid problems associated with the big 
bang singularity; in effect we are assuming that the scalar 



2m C\ m \ a+1 



(4.4) 



is a dimensionless parameter that determines the 
strength of the scalar field's contribution to the mass. 
We recall from Eqs. (|3.1|l and l|3.2() that a — —I when a 
is negative, and a = I + 1 when a is positive. We recall 
also that when a < 0, both rj and rjo are negative, which 
implies that (770/77) 2 > 1; when a > we have instead 
(Vo/v) 2 < 1- 

Because the hypergeometric function that appears in 
Eq. (|4.3|l is a simple polynomial of degree I — 1, closed- 
form expressions for the mass function can easily be ob- 
tained for small values of a. We list a few examples in 
Table II, and Figs. 1 and 2 provide graphs of the selected 
mass functions. 

The plots of Fig. 1 suggest that when a is a posi- 
tive integer, the mass function decreases at early times, 
but then increases until it is restored to its original 
value. Indeed, it is not difficult to show that the prop- 
erty m(rj)/mo — ► 1 as 77/770 — ► 00 follows directly from 
Eq. H4.3f) when a > 0. The plots of Fig. 2, on the other 
hand, suggest that when a < —2, only a fraction of the 
initial mass is eventually returned to the particle. This 
is confirmed by a similar calculation, which reveals that 
when a < 0, 



m(rj) 
mo 



1 - 



2 (21-1)!! 
Z — 1 (2Z — 2)!! 



(4.5) 



as 77/770 — > (recall that I = —a). When a = — 1 (de 
Sitter spacetime) the mass steadily decreases, and the 
limit goes to —00. 
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TABLE I: Reduced Green's function V(r], x; 77', x') for selected values of I, expressed in terms of R — \x — x'\. These expressions 
are obtained from Eq. I|3.10|l . 



I a V_ = 

{0,1} 

1 {-1,2} (wr 1 

2 {-2,3} | (r7 2 + r? ,2_ fl 2 )(W) - 2 

3 {-3,4} |[5r? 4 + 6 V 2 V ' 2 + 5t/ 4 - 10(t7 2 + r/ 2 )i? 2 + 5.R 4 ](W)" 3 

4 {-4,5} ^tf + v >i- R *)[ 7v i + 2r 1 2 r 1 ' 2 + 7 V ' 4 -U(r 1 2 +T 1 ' 2 )R 2 + 7R 4 ](r ] r 1 ')- i 




0.88 - 

0.86 ' 1 1 1 1 

1 2 3 4 5 6 

proper time: t/t 

FIG. 1: Mass function for selected positive values of a, plotted 
as a function of proper time t/tq = {v/rio) a+1 - The mass 
initially decreases but it returns to its initial value as r — > oo. 
The graph was obtained from Eq. I|4.3[l by setting £ = 0.1. 




1 2 3 4 5 6 

proper time: t/t 

FIG. 2: Mass function for selected negative values of a, plot- 
ted as a function of proper time t/tq. When a ^ — 1, 
t/to — (r//rio) a+1 , and the mass initially decreases but re- 
turns to a fraction of its initial value as r — > oo. When 
a = — 1 (de Sitter spacetime), t/tq = 1 — 111(77/770), and the 
mass decreases monotonically. The graph was obtained from 
Eq. JOJ by setting ( = 0.1. 



TABLE II: Mass function for selected values of a, expressed 
in terms of t = rj/r^Q. When a > t increases from 1 to 00. 
When a < t decreases from 1 to 0. These expressions are 
obtained from Eq. 14.31 . 



a 


m(rf)/mo = 1 — Q (• • •) 


2 


21n(t)t -3 


3 


3t" 3 - 3t~ 5 


4 


ft- 3 + l\n(t)t- 5 -ft- 7 


5 


35 .-3 , 45j.-5 45j.-7 35j.-9 
24 i t gl 8 l 24 ' 


-1 


-21n(t) 


-2 


3-3t 2 


-3 
-4 


f - f \n{t)t 2 - ft 4 

35 , 45 ,2 45,4 35 ,6 
24 + 8 1 8 1 24 L 



The classical theory of point scalar charges allows the 
particle's inertial mass to become negative. This may 
happen when the parameters a, q, ttlq, and 770 are suit- 
ably set (\a\ or C must be sufficiently large). A more com- 
plete theory would provide a description of what happens 
when m goes to zero. One possibility is that the parti- 
cle and its scalar charge simply cease to exist when the 
mass is completely radiated away; this scenario is pos- 
sible because the theory does not demand scalar-charge 
conservation. Another possibility is that the particle con- 
tinues to exist as a massless charged object [ll|. We shall 
remain neutral on the issue of continuing the life of the 
scalar charge beyond m = 0. 

V. SELF-FORCE ON A SLOWLY MOVING 
SCALAR CHARGE 

In this section we allow the scalar charge to move 
slowly relative to the cosmological fluid, we calculate the 
self-force that is acting on the particle, and we integrate 
the equations of motion to determine the particle's world 
line. To keep technicalities at a minimum we adopt a 
matter-dominated cosmology and choose a = 2, and we 
impose the condition \v\ <§; 1, where v — dz/drj is the 
particle's coordinate velocity vector. The metric of the 
spacetime is given by Eq. (|2.1|l with 0(77) = C77 2 . 

The slow-motion assumption implies that the particle's 
velocity vector can be approximated by u* 1 = (1, v)/a(rj). 
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The spatial components of the covariant acceleration are 
then 



Du a 
~dV 



1 



dv a 
drj 



V 



(5.1) 



The velocity of a particle moving on a geodesic there- 
fore satisfies dvjdr\ = —2v/r], and in time the particle is 
driven toward the same state of motion (v = 0) as the 
cosmological fluid. From a Newtonian point of view we 
might say that the cosmological expansion exerts a drag- 
ging force on the particle which compels it to come to 
an eventual rest. We shall adopt this point of view and 
refer to the 2v jr\ term in the equations of motion as the 
"cosmological dragging force" . 

The self-force acting on a scalar charge can be decom- 
posed into "Ricci" and "tail" forces. The Ricci force ap- 
pears as the first term on the right-hand side of Eq. I|2.9[) . 
while the tail force involves an integral over the past por- 
tion of the world line. For a matter-dominated cosmology 
the Ricci tensor is given by R x x = R y y = R z z = —R n „ = 
6/(C 2 r] 6 ), and the spatial components of the Ricci force 
are 



Ricci 



1 2q 2 

a 2 cy 



(5.2) 



Notice that the Ricci force is directed in the direction 
of v. For a matter-dominated cosmology the retarded 
Green's function is given by G = V/[a(r))a(r)')] with V = 
1/(7777'), so that G = C~ 2 (rir]')~ 3 . Evaluation of the tail 
integral on the right-hand side of Eq. Ij2.9|l gives 



rpa 

-'tail — 2 



392 i„JU. 



no 



(5.3) 



for the tail force. Notice that we have again truncated 
the tail integral so that it begins at 77 = 770. Notice also 
that the tail force is directed in the direction opposite to 
v. 

Combining Eqs. H5.1() - (|5.3() gives us the particle's equa- 
tions of motion, 



dv 
drj 



2 

—v 

V 



Crf"m(rj) 



31n^ 



(5.4) 



The first term on the right-hand side of Eq. I|5.4|) is the 
cosmological dragging force, which is directed along —v. 
The second term is the Ricci force, and this is directed 
along v. Finally, the third term is the tail force, and this 
is directed along —v. The sign of the self- force depends 
on time: The Ricci force dominates at early times (77/770 
close to unity) and it pushes on the particle, while the tail 
force dominates at late times and pulls on the particle. 

We shall assume that q 2 is small and that the self-force 
produces only a small correction to the particle's motion. 
In this regime we can replace 777,(77) by too on the right- 
hand side of Eq. (|5.4(1 . and our precise requirement is 
that 



q- 



2m C% 



« i; 



(5.5) 



notice that unlike q 2 , C is dimensionless. Introducing 
< = 77/770 > 1 we rewrite Eq. (|5.4|) as 



dv 2 
— = — v 
dt t 



^j(2-3hit)v. 



(5.6) 



The largest term on the right of Eq. (|5.6|1 is the cosmolog- 
ical dragging force, and this is only slightly offset by the 
self- force. Integration of Eq. (|5.6() . in the regime C "C 1> 
produces 



«(*) = 



t 2 



1 + S( <3 ~ 1 + 31ll< 



o(C 2 



(5.7) 



This result shows that when a charged (( 7^ 0) particle 
and an uncharged (£ = 0) particle start with the same 
initial velocity v(l), the charged particle moves faster 
at times t > 1. Ultimately, however, both particles are 
brought to rest by the cosmological expansion. 



VI. SELF-FORCE ON A SLOWLY MOVING 
ELECTRIC CHARGE 

That a charged particle moves faster than an un- 
charged particle in a cosmological spacetime (starting 
with equal initial velocities) is not an artifact of the 
scalar theory. The same conclusion applies in the con- 
text of electrodynamics. In this section we consider a 
point particle of (constant) mass to and electric charge e 
that moves slowly relative to the cosmological fluid. We 
once more calculate the self-force and integrate the par- 
ticle's equations of motion. Here we are able to select a 
more general class of cosmological models, and we let the 
scale factor take its general form 0(77) = Crj a with a no 
longer restricted to integers; we do, however, restrict a 
to be positive in order to keep the discussion simple. We 
once more impose the condition \v\ <C 1. 

It is known that the curved-spacetime version of 
Maxwell's theory is conformally invariant (see, for exam- 
ple, Appendix D of Ref. Q), and the metric of Eq. (f^TTf) 
is clearly conformally related to flat spacetime. This im- 
plies that the electromagnetic field produced by the mov- 
ing charge propagates without tails, and that the self- 
force is purely local. The charge's equations of motion 
were first determined by Hobbs [12j [. who corrected the 
earlier pioneering treatment of DeWitt and Brehme [Iflj 
(see also the reviews |2( and |3| for self-contained deriva- 
tions). In a conformally- flat spacetime, and in the ab- 
sence of an external force, the particle moves according 
to 



to-— = -e 2 (n + u^u^R^u* 
dr 3 



(6.1) 



here the self-force contains only a Ricci term, and the tail 
force is absent. 

The nonvanishing components of the Ricci tensor are 
R\ t = -ZaG- 2 T]- 2a - 2 and R x x = R\ = R z z = 
a(2a— l)C~ 2 rj~ 2a ~ 2 . In the slow-motion approximation 
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the particle's velocity vector is u M = (l,v)/a(rf), where 
v = dz/drj is the coordinate velocity. Using this infor- 
mation, we find that the spatial components of Eq. (|f> . 1|) 
reduce to 



dv a 2 .el 

~T = v + n a ( a + 1) ^ To 

d-q Tj 3 J mCri a + 2 



v. 



(6.2) 



The first term on the right-hand side of Eq. I|6.2[l is the 
cosmological dragging force; it is directed along — v. The 
second term is the electromagnetic self-force; it is di- 
rected along v and it pushes on the particle. 
In order to integrate Eq. I|6.2[l we rewrite it as 

dv a 4 , . > 1 ,„ . 

re- introducing t = 77/770 and defining the small dimen- 
sionless parameter 



c 



a+l ■ 



(6.4) 



Integration of Eq. (|6.3|) for £ <C 1 is straightforward, and 
we obtain 



v(t) 



v(l) 



+ ^r(* a * 1 -i)+°« a ) 



(6.5) 



This result leads to the same conclusion as in Sec. [V] 
When a charged (£ 7^ 0) particle and an uncharged 
(C = 0) particle start with the same initial velocity v(l), 
the charged particle moves faster at times t > 1. Both 
particles are eventually brought to rest by the cosmolog- 
ical expansion. 



VII. DISCUSSION: A PUSHING SELF-FORCE? 

It is somewhat surprising that a self-force (scalar or 
electromagnetic) would act so as to push the particle; 
this goes against a naive expectation that the self-force 
should always give rise to radiative damping. The exam- 
ples worked out in the preceding sections illustrate that 
this expectation is indeed naive. But there is no need to 
resort to curved-spacetime examples to show this. The 
expectation is naive also in the simplest context of slow- 
motion electrodynamics in flat spacetime. 

Consider a particle of mass m and electric charge e 
moving in the field of a charge — Q; we take e and Q 
to be positive quantities, and we suppose that Q is held 
fixed at the origin of a Cartesian coordinate system. The 
external force acting on e is the Coulomb force exerted 
by Q: 



-k—r, 



where k = (47re ) 1 , r 2 — x 2 + y 2 
distance from Q, and r = x/r is a radial unit vector; the 



(7.1) 



z 2 is the squared 



force of Eq. I|7.1|l is attractive. The self-force acting on e 
is 



-self 



Mo e • 

n r CXt 1 

D7TC m 



(7.2) 



where an overdot indicates differentiation with respect 
to time. This expression for the self-force is derived, for 
example, in Sec. 75 of the Landau and Lifshitz textbook 
fl4| . It is usually expressed in terms of the rate of change 
of the particle's acceleration vector, but Landau and Lif- 
shitz explain that it is consistent to make the replacement 
a — > i^ cx t /m on the right-hand side of Eq. I|7.2|l . A first- 
principles derivation of this equation is also presented in 
Ref. 

From Eq. I|7.1|l we calculate that the rate of change of 
the external force is given by .F cx t = —keQ(v — 3fr)/r 3 , 
where r = v • f is the radial component of the velocity 
vector. We decompose this vector in terms of radial and 
tangential components, 

v = v±+rf, v± ■ f = 0, (7.3) 

and we substitute this result into Eq. (|7.2|1 . This gives 



fi k e 3 Q , 
-Tsclf = (-v± 

OTTC 771 



2rf) 



(7.4) 



for the self-force acting on the moving charge. 

Suppose that the charge e is moving tangentially, so 
that r = 0. Then the self-force is directed along — v± = 
—v and we find that it is indeed a dragging force. Now 
suppose that the charge is moving radially inward, so that 
v± = and r < 0. Then the self- force is directed along 
— r, that is, in the same direction as the external force. 
In this situation the self-force pushes on the particle, in 
violation of the naive expectation. 

More generally, the rate of work done by the self-force, 
-F sc if • v, is proportional to — v 2 ^ + 2f 2 , where v± = \v±\. 
When v± > y/2r the work done on the particle is negative 
and the self-force exerts a drag. When, however, v± < 
\[2r the work done is positive and the self-force pushes. 

Does a pushing self-force violate energy conservation? 
The answer is in the negative, in spite of appearances 
to the contrary. Quite generally, the rate of work done 
on a particle of charge e by the self-force (which we now 
express in terms of the particle's acceleration a) is 



r fJ-oe 2 . fJ- e 2 

F se \{ ■ v = — — a ■ v = 



67TC 



67TC 



d 

dt {a - v) 



The last term on the right-hand side gives the electro- 
magnetic power radiated by the particle, and we have 
the statement 



[rate of work done] + [radiated power] 



67TC dt 



[a-v). 
(7.5) 

Integrating this equation over a finite time interval gives 
rise to the following conservation statement: 

[work done] + [radiated energy] = 

[total change in ^e 2 {a ■ v)/(6ttc)]. (7.6) 
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The term on the right can be interpreted as a change in 
the system's "internal energy" as it moves from one con- 
figuration to another. (The system is the charged particle 
plus the agent responsible for the external force.) This 
form of energy is associated with the part of the electro- 
magnetic field that stays bound to the system instead of 
propagating freely as electromagnetic radiation. 

In some circumstances the total change in internal en- 
ergy is zero. For example, when the motion is peri- 
odic the total change vanishes after one complete pe- 
riod. As another example, the total change vanishes 
when the motion proceeds from a state of no accel- 
eration to another state of no acceleration. In such 
circumstances the statement of Eq. ()7.6|) simplifies to 



[work done] + [radiated energy] = 0, and the work done 
by the self-force must be negative. In more general cir- 
cumstances, however, the statement of energy conserva- 
tion must include the change in internal energy, and it 
becomes possible for the self-force to do positive work 
on the particle. This was well illustrated in the previous 
discussion. 
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